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Abstract. Relativistic field equations for a gas in special and general relativity are determined from the Boltzmann equation. 
The constitutive equations are obtained from the Chapman-Enskog methodology applied to a relativistic model equation 
proposed by Anderson and Witting. Two applications in general relativity are considered: one refers to a gas in a homogeneous 
and isotropic Universe where irreversible processes are present during its evolution; in the other it is analyzed a gas under the 
influence of a spherically symmetrical non-rotating and uncharged source of the gravitational field. 



Keywords: Boltzmann equation, special relativity, general relativity 
PACS: 51.10.+y, 05.20.Dd, 03.30.+p, 04.20.-q 



INTRODUCTION 

We may say that the development of the relativistic kinetic theory goes back to 1911 when Jiittner [1] succeeded to 
derive an equilibrium distribution function for a relativistic gas, which in the non-relativistic limiting case reduces to 
the Maxwellian distribution function. The relativistic Bose-Einstein and Fermi-Dirac distribution functions were also 
derived by him in 1928 [2]. At the beginning of the 1940s the covariant formulation of the Boltzmann equation was 
proposed by Lichnerowicz and Marrot [3], while the determination of the transport coefficients from the Boltzmann 
equation by using the Chapman-Enskog method was obtained by Israel [4] and Kelly [5] in the 1960s. The general 
relativistic formulation of the Boltzmann equation goes back also to the sixties of the last century when Chernikov [6] 
analyzed a relativistic gas in a gravitational field. 

The aim of this work is to analyze the Boltzmann equation within the framework of special and general relativity 
and to derive the corresponding macroscopic field equations. 

First the Boltzmann equation in special relativity is analyzed and the balance equations for the particle four-flow 
and energy-momentum tensor are obtained. As in the case of non-relativistic fluid dynamics these equations are useful 
to describe a relativistic fluid with five scalar fields of particle number density, four-velocity and temperature. The 
constitutive equations for the non-equilibrium pressure, pressure deviator and heat flux, as well as the corresponding 
transport coefficients of bulk and shear viscosities and thermal conductivity, are determined by using the Chapman- 
Enskog methodology applied to the Anderson and Witting model of the Boltzmann equation [7]. From the knowledge 
of the constitutive equations, the system of balance equations turns out to become a system of field equations for the 
basic fields and problems related to relativistic fluids can be analyzed. 

Next, relativistic fluids in general relativity are analyzed within the framework of Boltzmann equation. In the first 
application a gas in a homogeneous and isotropic Universe is considered where irreversible processes are present 
during its evolution. The irreversible processes are characterized by a non-equilibrium pressure which is determined 
from the Anderson and Witting model by using the Chapman-Enskog methodology. Another application refers to a gas 
- which is not the source of the gravitational field - under the influence of a Schwarzschild metric, which is the solution 
of Einstein's field equations for a spherically symmetrical non-rotating and uncharged source of the gravitational field. 
For this case, the expression for the momentum density balance equation is obtained in the non-relativistic limiting 
case and for weak gravitational fields. 

In this work Greek indices run from a = 0, 1 , 2, 3, while Latin indices from i = 1,2,3. 



BOLTZMANN EQUATION IN SPECIAL RELATIVITY 



We consider a single non-degenerate relativistic gas in a Minkowski space with metric tensor r\ a p = 
diag (1, — 1, — 1, — 1), where a particle of rest mass denoted by m is characterized by the space-time coordinates 



(x a ) = (x = ct,x) and by the momentum four-vector (p a ) = (p ,p), where c is the speed of light. The momentum 
four- vector has a constant length mc which implies that p° can be expressed in terms of p by p° = \/IpI 2 +m 2 c 2 . 

The state of the relativistic gas - in the phase space spanned by the space-time and momentum coordinates - is 
characterized by the one-particle distribution function f(x a ,p a ) = /(x,p,f), where the quantity /(x,p,f)rf 3 jcJ 3 p = 
/ (x 7 p,t)dx l dx 2 dx 3 d p l d p 2 d /? 3 gives the number of particles at time t in the volume element d 3 x about x and with 
momenta in a range c/ 3 p about p. 

In relativity we have to take into account the behavior of the quantities under Lorentz transformations. The volume 
element dr'xcPp is an invariant under Lorentz transformations, i.e., dr'xcPp = d^x'd 3 p'. Furthermore, the number of 
particles in a volume element is a scalar invariant, since all observers will count the same number of particles. Hence, 
we may conclude that the one-particle distribution function f(x, p, t ) is also a scalar invariant. It can be proved that the 
volume element d 3 p is not an invariant but only the ratio d? p' / p' = d 3 p/po- 

The evolution of the one-particle distribution function in the phase space is ruled by the Boltzmann equation, which 
in the absence of external forces reads 

In the above equation we have introduced the abbreviations fl = /(x,p*,£), /' = /(x,p',f), /» = /(x,p„,,?), / = 
/(x,p,f)„ where p and p» denote the momenta of two particles before a binary collision and p' and p'„ are the 
corresponding momenta after collision. The pre and post collisional momentum four-vectors are connected by the 
energy-momentum conservation law p a + p" = p' a + p'". Furthermore, F = \J (p"pa) 2 — m 4 c 4 is the invariant flux, 
which in the non-relativistic limiting case it is proportional to the modulus of the relative velocity. The differential 
cross-section and the element of solid angle that characterize a binary collision are denoted by a and dD., respectively. 

The multiplication of the Boltzmann equation (1) by an arbitrary function y(pP ) and integration of the resulting 
equation with respect to the invariant d 3 p/ po, yields the following general equation of transfer 

A / W a f— = J f (Y+ Y*-V- <) if J' ~ Uf)FadQ^^-, (2) 
dx a J po 4 J p*o po 

which is only a function of the space-time coordinates. As in the non-relativistic case the right-hand side of the above 
equation is obtained by evoking the symmetry properties of the collision operator. When the underlined term in (2) 
vanishes, l// becomes a summational invariant. The necessary and sufficient condition to be a summation invariant is 
that \j/ must be given by y/ = A + B a p a , where A is an arbitrary scalar and B a an arbitrary four- vector that do not 
depend on the momentum four- vector p a . 

The macroscopic description of a relativistic gas may be characterized by the two first moments of the distribution 
function, namely 

N a = c j p a f—^-, T<*P= c r p "pPf d -^, (3) 
J po J po 

which are known as the particle four-flow and the energy-momentum tensor, respectively. Their balance equations are 
obtained from the transfer equation (2) by choosing y/ = c and \jf = cpP, yielding 

d a N a =0, d a T a P=0, (4) 

respectively. They represent the conservation laws of the particle four-flow and of the energy-momentum tensor. 

In the analysis of relativistic fluids it is usual to decompose the particle four-flow and the energy-momentum tensor 
in terms of quantities that appear in non-relativistic fluid dynamics, like particle number density «, energy per particle 
e, hydrostatic pressure p, non-equilibrium pressure GJ, pressure deviator P"^ (the traceless part of the pressure tensor) 
and heat flux q a . There exist two decompositions for the particle four-flow and energy-momentum tensor as functions 
of the above mentioned quantities and they are known as the decompositions of Eckart and of Landau-Lifshitz (see 
e.g. [8]). Here we shall adopt the Landau-Lifshitz decomposition, namely 

N a = nU a -^-, T a P =F a P -(p + ®)A a P + ^U a U^. (5) 
h c z 

The above decomposition is written in terms of the four-velocity U a (such that U a U a = c 2 ) and of the projector 
A aj g = ?7 a j3 — U a Up/c 2 (such that A^Up = 0). Furthermore, we have introduced the enthalpy per particle h = e+p/n. 



As in the non-relativistic case, the equilibrium is characterized by the condition /£/' = f*f which renders a 
vanishing collision term of the Boltzmann equation (1). This last relationship implies that In/* + ln/* = In/* + ln/ 
and In/ represents a summational invariant whose expression is given by In/ = A +B a p a - The terms A and B a are 
determined from the equilibrium values of the particle four-flow and energy-momentum tensor and it follows that the 
equilibrium distribution function is represented by the Maxwell- Juttner distribution function 

denoting the modified Bessel function of second kind. Note that the distribution function depends on £ = mc 2 /kT 
which represents the ratio of the rest energy mc 2 of a particle and the thermal energy of the gas IcT, where k is the 
Boltzmann constant and T the temperature of the gas. The limiting case £ ^> 1 represents a non-relativistic regime of 
the gas, while £ <C 1 its ultra-relativistic regime. For example, for an electron of rest mass m w 9.1 x 10~ 31 kg we 
have that £ = 5.9 x 10 9 /T, indicating that only for very high temperatures the electron might be considered in the 
relativistic regime. 

If we are interested in the description of the relativistic gas as a five-field scalar theory characterized by the fields 
of particle number density «, four- velocity U a and temperature T, we have to determine the constitutive equations for 
the energy per particle e, hydrostatic pressure p, non-equilibrium pressure GS, pressure deviator P"^ and heat flux q a 
as functions of (n,U a ,T). The determination of the constitutive equations for the energy per particle and hydrostatic 
pressure follows from the insertion of the Maxwell-Juttner distribution function (6) into the definition of the energy- 
momentum tensor (3)2 and the integration of the resulting equation, yielding 



mc 2 



k 2 (Q C 



p = nkT. (7) 



Hence, the expression for the hydrostatic pressure is the same as the one for a non-relativistic gas. In the non-relativistic 
and ultra-relativistic limiting cases the energy per particle becomes 



mc 2 + — kT 



3. 



e = 3kT + . . . , (8) 



respectively. Note that the first term in the non-relativistic limiting case (8)i is the rest energy of a particle mc , the 
second one 3kT /2 represents the classical value of the energy density, while the third one is a relativistic correction. By 
considering only the first two terms in this expansion it follows the well-known relationship between the pressure and 
the energy density of a non-relativistic gas, namely, p = 2n(e — mc 2 )/3. Furthermore, in the ultra-relativistic limiting 
case the following relationship p =ne/3 holds. 



Model equation and Chapman-Enskog method 

The constitutive equations for the non-equilibrium pressure GJ, pressure deviator P"^ and heat flux q a can be 
obtained from the Boltzmann equation by using the Chapman-Enskog methodology. Here we shall use the simple 
model equation proposed by Anderson and Witting [7], which is a generalization of the BGK model for a relativistic 
gas. The model equation for the Boltzmann equation in this case is written as 

p a (U a Df + W a f) = - l ^(f-f { % (9) 

where T is a characteristic time of order of the mean free time. Above we have introduced the differential operators 
D and V a through the decomposition d a = (U a D/c 2 + V«). The differential operators D and V a reduce to a partial 
time derivative and a spatial gradient in a Lorentz rest frame, respectively. 

For processes close to equilibrium we may write / = /W (1 +0), where the deviation <j) from the Maxwell-Juttner 
distribution function is considered as a small quantity. Now by using the Chapman-Enskog methodology we insert the 
equilibrium distribution function on the left-hand side of (9) and the representation / = /'°' (1 + <j>) into its right-hand 
side and get 

p«(c/ a D/( o )+V a /(°))=-^0. (10) 



Next we perform the derivatives which occur in the left-hand side of above equation and eliminate the derivatives Dn, 
DT and DU a by using the balance equations for the particle number density, energy density and momentum density 
of an Eulerian relativistic gas (where G7 = 0, P"^ = and q a = 0), namely, 



Dn + nV a U a = 0, nc v DT + pV a U a = 0, 



n -^DU a = V V 



(11) 



Hence, the deviation from the Maxwell-Juttner distribution function can be determined, yielding 
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(12) 



In the above equation we have introduced the heat capacity per particle at constant volume c v and the traceless part of 
the velocity gradient W a UP defined by 



AJAJ 
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y<5 



(13) 



Note that the deviation is a function of the thermodynamic forces which are underlined in (12): the divergence of 

velocity V a U a , the traceless part of the velocity gradient V a ljP and a combination of the temperature and pressure 
gradients (y a T — ^J a p)- The contribution of the pressure gradient refers only to a relativistic correction. 

The constitutive equations for the non-equilibrium pressure, pressure deviator and heat flux can be obtained now 
from the insertion of / = /'°'(1 + 0) together with (12) into the definitions of the particle four-flow and energy- 
momentum tensor (3) and integration of the resulting equations. Hence it follows 



m = -r]V a U a , P a/3 =2 J uV a t// 3 , q a = X 



V a T-—V a p 
nh 



(14) 



The above equations represent the laws of Navier-Stokes and Fourier and we may identify the scalar coefficients 77, 
and A as the coefficients of bulk viscosity, shear viscosity and thermal conductivity, respectively. Their expressions 
read 
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(15) 
(16) 



where Kii denotes the following integral for the modified Bessel function Kii (C) = Io° e ^ cosh ' dt / cosht . 
In the non-relativistic limiting case (£ ^> 1 ) the transport coefficients become 
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while in the ultra-relativistic limiting case they read 
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(17) 



(18) 



In the non-relativistic limiting case the ratio A/jU ~ 5k /2m has the same value as that of the BGK model. Furthermore, 
we note that the bulk viscosity coefficient is important only in the relativistic case, since it vanishes in the non- 
relativistic and in the ultra-relativistic limiting cases. 

Once the constitutive equations (7) and (14) are known, the balance equations (4) become a system of field equations 
for the five scalar fields of particle number density n, four-velocity U a and temperature T and problems concerning 
relativistic fluids in special relativity could be analyzed. 



BOLTZMANN EQUATION IN GRAVITATIONAL FIELDS 



For the analysis of a relativistic gas in the presence of a gravitational field we have to consider a Riemannian 
space characterized by the metric tensor g^. Here the invariant elements in the phase space are \/—g 1 d ? 'p'/p' Q = 
yf—gcPp I po and gifixcfi p /po = g'cPx'cPp' j p' , where g = det(^ jUV ) is the determinant of the metric tensor. Moreover, 
the constant length of the momentum four- vector g^vP^p v = m 2 c 2 implies that 

P° = Po - gOiP'/ goo , with po = g QO m 2 c 2 + (g 0i g Q ] - g QO gjj ) p'pi . (19) 
The Boltzmann equation in gravitational fields is written as 

- = J M -f.f)F°^^, (20) 

where F° v = g az (d^g^ + d v g^ x — dxg^v) /2 are Christoffel symbols. 

Following the same methodology of the special relativistic case, we may define the two first moments of the 
distribution function and obtain their respective balance equations from the Boltzmann equation (20), which read 

N^= I >/V=S — , A^=<9 M M<+rpV CT =0, (21) 
J po 

T^ v = I P^p V fV^8—, ^ v ;M = a M T^ v +T%T av +Tl a n° =0. (22) 
" Po 

In the above equations the semicolon refers to a covariant derivative. 

For the solution of a problem concerning a relativistic gas in the presence of a gravitational field, it is necessary to 
determine the ten components of the metric tensor g^ v , which is achieved by using Einstein's field equations 

*Mv-|«gMv = -^r MVj where fy v = dyT^ - d T r£ v + r^r^ - rjtyl^, , R = g» v R liv , (23) 

represent the Ricci tensor and the curvature scalar, respectively. Moreover, G denotes the gravitational constant. In 
the right-hand side of Einstein's field equations the energy-momentum tensor represents the sources of the gravi- 
tational field, while its left-hand side corresponds to the geometry of the curved space-time. The conservation of 
the energy-momentum tensor T^.^ = follows also from Einstein's field equation thanks to the Bianchi identity 



Relativistic gas in a Robertson- Walker metric 

The search for cosmological solutions of our Universe are founded on Einstein's field equations and on the 
hypothesis that the Universe is spatially homogeneous and isotropic at large scales, i.e., of order of 100 Mpc w 
3.1 x 10 24 m w 3.26 x 10 8 light-years. The isotropy and homogeneity assumptions imply that the line element ds 
is given by the Friedmann-Lamaitre-Robertson- Walker (FLRW) metric, which for a spatially flat Universe reads 

ds 2 = (c 2 dt} 2 — a(t) 2 (dx) 2 + (dy) 2 + (dz) 2 ■ This metric is a function of only one unknown parameter which is 
the cosmic scale factor a(t). In addition from the hypothesis of homogeneity and isotropy of the space-time, the 
pressure deviator P ^ and of the heat flux q a must be absent from the energy-momentum tensor, and in this case it 
reduces to T^ v = diag(e«, — p — CJ, —p — tfi,—p — d). 

For the FLRW metric the non-vanishing components of the metric tensor and of the Christoffel symbols are given 
by goo = 1, gjj = ~a 2 8ij, g = —a 6 , = dadjj, T' j = d'jd/a, where the dot denotes the derivative with respect to 

x° = ct . 

The components 00 and 1 l(say) of the Einstein's field equations (23) i lead to two equations which may be used to 
analyze cosmological problems. These equations - called the Friedmann and the acceleration equations - read 

d\ 2 8ttG a \%G 

- =TX e "> -=-TT [en + 3(p + ® )}, (24) 
a I 3c 4 a 3c 4 



respectively. In order to obtain a closed system of equations for the cosmic scale factor a and energy density en 
it is necessary to close the system of equations (24) with constitutive laws for the hydrostatic pressure p and for 
the non-equilibrium pressure UJ. The thermal equation of state p = p(en) gives the dependence of the pressure on the 
energy density. To obtain the constitutive equation for the non-equilibrium pressure, we shall use the Chapman-Enskog 
methodology for the Anderson and Witting model (20). In the case of a homogeneous and isotropic Universe we may 
neglect the spatial gradients, so that (20) can be written as 
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(25) 



by considering a spatially flat Universe and a comoving frame where (U ll ) = (c,0). In the above equation the Maxwell- 
Juttner distribution function is given by 



/ (0) (x,P,0 = 
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kT 



(26) 



Note that f(°> depends on the cosmic scale factor a. 

By following the same methodology above, we write / = /'°' ( 1 + ) and insert this representation on the right-hand 
side of (25) and on its left-hand side we perform the derivatives of /W. Hence, we get 



= -CT 
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(27) 



The second equality above was obtained through the elimination of the time derivatives of the particle number density 
and of the temperature by considering the equations 



h + 3n- 
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3kT d 
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0, 



(28) 



which follows from (21)2 and from the Oth-component of (22)2 with CT = 0, respectively. 

The following projection of the energy-momentum tensor corresponds to the sum of the hydrostatic pressure with 
the non-equilibrium pressure 
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(29) 



In order to perform the above integral we introduce a new variable y = po/mc = \J 1 +a 2 |p| 2 / (mc) 2 so that we can 
write 

i i , . mc ydy 



and 



dp 



« (j 2 -l)2 

Next, the insertion of f — /^(l + <j>) together with (27) into (29) and by taking into account (30), we get 



-GJ = 
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(30) 



(31) 



From the integration of the above equation we get the following constitutive equation for the non-equilibrium pressure: 



US = -3TJ- 



(32) 



We identify 77 in (32) as the coefficient of bulk viscosity, since its expression is the same as that given by (15). 
Furthermore, from the comparison of (32) with the constitutive equation for the non-equilibrium pressure - given in 
terms of the divergence of the four- velocity ( 14) 1 - we conclude that 3d /a plays the same role as the divergence of the 
four-velocity. Owing to the fact that the bulk viscosity is a positive quantity the non-equilibrium pressure decreases 
when the Universe is expanding (d > 0) and it increases when it is contracting (d < 0). 

Now a cosmological problem can be solved from the system of coupled differential equations (24) together with 
the constitutive equations for the pressure p — p(en) and for the non-equilibrium pressure (32), since in this case (24) 
becomes a system of field equations for the determination of the cosmic scale factor a(t ) and energy density ne{t ). 



Relativistic gas in a Schwarzschild metric 

Now we shall analyze a gas under the influence of a Schwarzschild metric, which is the solution of Einstein's field 
equation for a spherically symmetrical non-rotating and uncharged source of the gravitational field. In terms of the 
spherical coordinates (f, 0,(p) the Schwarzschild metric reads 



ds l = 1 - 



2GM\ 



dr 2 - f 2 ( dB L + sin z 9d(p' 



(33) 



where M is the total mass of the spherical source. If we introduce a new radial coordinate r = djjx'x-i through the 

. 2 

GM 



(34) 



relationship f = r y \ + , we get the isotropic Schwarzschild metric in the new radial coordinate 
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For the isotropic metric (34), the Christoffel symbols read 

r° = o, r?. = o, ^ = 0, rf. = o (i^j^k), r°=-i 
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where the underlined indices are not summed. Furthermore, it follows that po = ^fgo~yJm 2 c 2 +g\\\)\ 2 , p° = Po/go, 



^/—g = J gag], while the components of the four- velocity in a comoving frame reads = (c/ ^/gj),0). 

Now by writing / = /'°'(l + 0) and following the Chapman-Enskog methodology we get from the Boltzmann 
equation (20) 
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The multiplication of (37) by ^J—gd^pj po and integration of the resulting equation leads to the balance equation of 
particle number density 

U v ^-+nU v - v =0. (38) 
dx v 

Furthermore, the multiplication of (37) by p^ y^gd^p / po and subsequent integration imply into the balance equations 
for the energy density and momentum density of an Eulerian fluid. They read 



dT 

nc v U v -=— - +pU v v = 0, — th component, 
dx v 
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(39) 
(40) 



respectively. 

It is interesting to discuss the balance equation for the momentum density (40), since the underlined term refers to 
the presence of the gravitational field. By introducing the gravitational potential 



GM 



d<t> GM xi 
so that 3-r = —5-6*; — , 
ax K r L r 



(41) 



we can write the balance equation for the momentum density (40) as 



nh jjvjjk d P vk, P *k) d * 



_8_ 2X3 «3 C 
C 2 + ZK 2 + K 2 \ 



= 0. (42) 



The ratio <t>/c 2 can be estimate at the surface of some bodies: 

1. Earth: M ffi 5.97 x 10 24 kg; # e 6.38 x 10 6 m; 4>/c 2 7 x 1(T 10 ; 

2. Sun: M Q w 1.99 X 10 30 kg; # w 6.96 x 10 8 m; 4>/c 2 w 2.2 x 1(T 6 ; 

3. White dwarf: M w 1.02M Q ; fi«5.4x 10 6 m; 4>/c 2 « 2.8 x 1(T 4 ; 

4. Neutron star: M f« M Q ; 7? w 2 x 10 4 m; 4>/c 2 7.5 x 1(T 2 ; 

5. Black hole: M 3M Q ; 7? 3 x 10 3 m; 4>/c 2 w 1.5. 

From the above estimates we infer that in most cases we can use the approximation <t>/c 2 <C 1. 

Hence, when <t>/c 2 <C 1 and by considering the non-relativistic limiting case £ 3> 1, the balance equation of the 
momentum density (42) becomes 





!_ 2 * + ™ +. ..! + _+-+... 5*^=0. (43) 
c 2 4c 4 / \ 2C 8C 2 / dxJ 



By neglecting the underlined terms the above equation reduces to the usual Newton's second law in the case of a 
non-relativistic gas in the presence of a weak gravitational field. The underlined terms represent the relativistic and the 
gravitational potential corrections. 



CONCLUSIONS 

In this work we have derived the field equations for relativistic gases in special and general relativity from the 
Boltzmann equation. The transport coefficients of a relativistic gas in special relativity were obtained from the 
Anderson and Witting model equation. In general relativity we have analyzed first a homogeneous, isotropic and flat 
Universe, where irreversible processes are taken into account and derived the corresponding field equations. Next, a 
relativistic gas under the influence of a spherically symmetrical non-rotating and uncharged source of the gravitational 
field was considered and the corresponding balance equations were obtained. 
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